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Abstract 

We review several of the normal state properties of the cuprates in an attempt to 
establish an organizing principle from which pseudogap phenomena, broad spectral 
features, T— linear resistivity, and spectral weight transfer emerge. We first show 
that standard field theories with a single critical length scale cannot capture the 
T— linear resistivity as long as the charge carriers are critical. What seems to be 
missing is an additional length scale, which may or may not be critical. Second, we 
prove a generalised version of Luttinger's theorem for a Mott insulator. Namely, 
regardless of the spatial dimension, the Fermi surface of the non-interacting system 
is converted into a surface of zeros of the single-particle Green function when the 
Mott insulator posesses particle-hole symmetry. Only in the presence of particle-hole 
symmetry does the volume of the surface of zeros equal the particle density. The 
surface of zeros persists at finite doping and is shown to provide a framework from 
which pseudogaps, broad spectral features, spectral weight transfer on the Mott gap 
scale can be understood. 



1 Introduction 



The normal state of the high-temperature copper-oxide superconductors (cuprates 
for short) is anything but normal First, all parent cuprates are Mott insulators pLj- 
Such insulators do not insulate because the band is full, but rather because 
strong local electron correlations dynamically generate a charge gap by splin- 
tering a half-filled band into lower and upper Hubbard bands. Second, once 
doped, Mott insulators still insulate. In fact, numerous experiments [2][3]fill5] in 
which a large magnetic field is applied to kill superconductivity reveal that the 
resistivity diverges as the temperature tends to zero below a critical doping 
level, Xc- Above Xc, traditional Fermi liquid behaviour is observed in the over- 
doped regime. Third, electron quasiparticles are nowhere to be found in the 
underdoped regime where the ubiquitous pseudogap [71l8|9] phenomenon ob- 
tains. Fourth, near optimal doping, the resistivity scales as a linear[TO] function 
of temperature as opposed to the quadratic behaviour indicative of a Fermi 
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liquid. Such oddities did not, of course, beset a theoretical understanding of 
superconductivity in metals such as aluminum because the normal state is 
accurately rooted in Fermi liquid theory. In the cuprates, the central problem 
of the normal state is understanding how the disparate phenomena just de- 
scribed arise from strong electron correlations intrinsic to the Mott state. Sim- 
ply stated, the question is: is there a general organizing principle that captures 
Mottness? Mottness refers to those features of doped Mott insulators that are 
not deducible from ordering. While it might be that ordering is deducible 
from Mottness, Mottness is distinct nomologically from ordering. Indeed, it 
is Mottness that makes the cuprate problem largely intractable. It is for this 
reason that numerous theories of the cuprates get Mottness clearly wrong, ig- 
nore it, or deny its existence. In this paper, we prove two theorems which bare 
directly on the nature of the low-energy theory for a doped Mott insulator. 
Both theorems point to physics beyond that which is captured by standard 
field theories based on one-parameter correlation length scaling. In the first 
theorem, we show explicitly that within single parameter scaling, quantum 
criticality cannot account for the linear T resistivity. Second, we show that for 
Mott insulators possessing particle-hole symmetry, the single-particle Green 
function vanishes at the Fermi surface of the non-interacting system. While 
the surface of zeros is shown to persist in the absence of particle-hole symme- 
try, only in the presence of such a symmetry does the volume of the surface 
of zeros equal the particle density. The surface of zeros provides a framework 
from which pseudogaps, broad spectral features, and spectral weight transfer 
on the Mott gap scale can be understood. 



2 T-Linear Resistivity 

One of the leading explanations for the T— linear resistivity is quantum criticality [TT][T^ . 
We do not analyse here specific proposals such as the highly successful marginal- 
Fermi liquid theory [13]. which is in principle consistent with a quantum crit- 
ical point, but rather dissect the assumptions of quantum criticality which 
lead to T— linear resistivity. At the quantum critical coupling, or in the quan- 
tum critical regime, ksT is the only energy scale governing collisions between 
quasiparticle excitations of the order parameter. Consequently, the transport 
relaxation rate. 



scales as a linear function of temperature, thereby implying a T— linear re- 
sistivity if (naively) the scattering rate is what solely dictates the transport 
coefficients. That temperature is the only scale in the quantum critical regime 
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Fig. 1. Heuristic phase diagram of the cuprate superconductors as a function of 
temperature and hole doping level. AF represents antiferromagnet, SG, the spin 
glass, SC the superconductor, PG the pseudogap regime in which the single-particle 
specrum develops a dip and FL the Fermi liquid. The spin-glass phase ceases at a 
critical doping level (quantum critical point, QCP) inside the dome. The dashed lines 
represent crossovers not critical behaviour. The strange-metal behaviour, T— linear 
resistivity, in the funnel-shaped regime has been attributed to quantum critical 
behaviour. A scaling analysis of the conductivity at the quantum critical point 
rules out this scenario, however. 

regardless of the nature of the quasiparticle interactions is a consequence of 
universality. Eq. ([1]) holds as long as the inequalities T > |A| and t < 1/|A| 
are maintained, A the energy scale measuring the distance from the critical 
point and t the observation time. The energy scale A oc S^" varies as a function 
of the tuning parameter 6 = g — Qc, where u is the correlation length exponent 
and z is the dynamical exponent. At the critical coupling 6 = or g = Qc, A 
vanishes. Ultimately, the observation time constraint, t < 1/|A| implies that 
only at the quantum critical point does the T— linear scattering rate obtain 
for all times. That the quantum critical regime is funnel-shaped follows from 
the inequality T > |A|. The funnel-shaped critical region should be bounded 
by a temperature Tuppor above which the system is controlled by high-energy 
processes. Whether or not quantum criticality is operative up to temperatures 
of order T = lOOO/T in the cuprates is questionable. 

We have recently formulated [H] a general scaling argument for the conductiv- 
ity to understand if T— linear resistivity is consistent with the general claims 
of quantum criticality. The argument is based on three general assumptions: 1) 
the critical degrees of freedom carry the current, 2) only a single length scale 
is relevant near the critical point and 3) charge is conserved. Granted true, 
these assumptions rule out quantum criticality as the cause of the T— linear 
resistivity. 

Before we present the proof, consider as a precursor to the Drude formula for 
the conductivity. 
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with (jjpi the plasma frequency and Tti the transport relaxation time. If we sub- 
stitute Eq. ([1]) into Eq. ([2]), T— linear resistivity follows in the zero-frequency 
limit. However, explicitly appearing in the formula for the resultant conductiv- 
ity is the Upi. While the plasma frequency is commonplace in typical transport 
theories of metals, we show below that such a degree of freedom does not exist 
in quantum critical theories. Consequently, the predictions of Drude theory are 
not compatible with those of quantum criticality. That inconsistencies lie with 
the Drude formula and quantum criticality can be seen immediately from the 
work of van der Marel[15j and co-workers who used Eq. ([2]) to collapse their op- 
tical conductivity data as a function oiuj/T. True quantum critical behaviour 
would be described by an optical conductivity of the form, T'^ fioj /T). Con- 
trastly, they find|31j that /i = 1 for uj/T < 1.5 and ~ 0.5 for uo/T > 3 if 
f{ijj/T) is described by the Drude formula. 

Consider a general action S whose microscopic details are unimportant as 
long as the current is critical. An externally applied electromagnetic vector 
potential A^, = 0, 1, . . . , rf, couples to the electrical current, j^, such that 

S ^ S + j dT(fxA^' j^. (3) 

Under the one-parameter scaling hypothesis for quantum systems, the spatial 
correlations in a volume smaller than the correlation volume, and temporal 
correlations on a time scale shorter than oc are small, and space-time 
regions of size ^'^^t behave as independent blocks. Using this hypothesis, we 
write the scaling form for the singular part of the logarithm of the partition 
function by counting the number of correlated volumes in the whole system: 

lnZ = ^F(5e'^^{AU'^^}), (4) 

In this expression, L is the system size, jS = l/fc^T the inverse temperature, 
6 the distance from the critical point, and ds and dA the scaling dimensions 
of the critical coupling and vector potential, respectively. The variables A\ = 
A''{uj = A^^r ) correspond to the (uniform, k = 0) electromagnetic vector 
potential at the scaled frequency A = ujC,t, and i = 1, . . . ,d labels the spatial 
components. Two derivatives of the logarithm of the partition function with 
respect to the electromagnetic gauge A^u), 
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determine the conductivity for carriers with charge Q. We have exphcitly set 
5 = as our focus is the quantum critical regime. At finite temperature, the 
time correlation length is cutoff by the temperature as 1/T, and . 
The engineering dimension of the electromagnetic gauge is unity (d^ = 1). 
From charge conservation, the current operators cannot acquire an anomalous 
dimension; hence, ci^ = 1 is exact |16j. We then arrive at the general scaling 
form 



for the conductivity where E is an explicit function only of the ratio, uj/T. 
(The ij tensor indices have been dropped for simplicity.) This scaling form 
generalizes to finite T and uj the T = frequency dependent critical conduc- 
tivity originally obtained by Wen [16]. The generic scaling form, Eq. ([6]), is 
also in agreement with that proposed by Damle and Sachdev[17j in their ex- 
tensive study of collision-dominated transport near a quantum critical point 
(see also the scaling analysis in Ref. [H]). What the current derivation lays 
plain is that regardless of the underlying statistics or microscopic details of 
the Hamiltonian, be it bosonic (as in the work of Damle and Sachdev[l7]) or 
otherwise, be it disordered or not, the general scaling form of the conductivity 
is unchanged. The Anderson metal-insulator transition md = 2 -(- e, which can 
be thought of as a quantum phase transition where the dimensionless disorder 
strength is the control parameter [201IT9] . obeys the scaling function derived 
here for the conductivity. 

At zero frequency, the dc limit, we obtain 



In general S(0) 7^ 0. Else, the conductivity is determined entirely by the 
non-singular and hence non-critical part of the free energy. The cuprates are 
anisotropic 3-dimensional systems. Hence, the relevant dimension for the crit- 
ical modes is = 3 not d = 2. In the latter case, the temperature prefactor is 
constant. For d = 3, we find that T— linear resistivity obtains only if z = — 1. 
Such a negative value of z is unphysical in standard commutative |24j field the- 
ories as it implies that energy scales diverge for long wavelength fluctuations 
at the critical point. 

While it is certainly reasonable to assume that S(0) is finite at zero temper- 
ature, it is certainly a possibility that S(0) might in fact diverge. A possible 
scenario of how this state of affairs might obtain is that a dangerously irrel- 
evant operator could govern the conductivity. In this case, S(0) ~ l/T^, and 
T— linear behaviour obtains if p = {d — 2)/z+l. This can only occur above the 
upper critical dimension. In this regime, all criticality is mean- field-like. Hence, 
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a possibility that cannot be eliminated at this time is that all criticality in the 
cuprates is inherently mean-field and dangerously irrelevant operators control 
the conductivity. For strongly correlated systems, however, it is unlikely that 
criticality is inherently mean-field. 

Consequently, three options are available at this point: 1) quantum criticality 
has nothing to do with the problem, 2) the current is carried by non-critical 
degrees of freedom or 3) new quantum critical scenarios in which additional 
length scales describe the physics. In a scenario involving non-critical degrees 
of freedom, fermionic charge carriers in the normal state of the cuprates could 
couple to a critical bosonic mode. Such an account is similar to that in mag- 
netic systems [11 J in which fermions scatter off massless bosonic density or spin 
fluctuations and lead to an array of algebraic forms for the resistivity [22ll50] 
ranging from T^/^ to T^/^ in antiferromagnetic and ferromagnetic systems, 
respectively. While disorder can alter the exponent |50]. T— linear resistivity 
results only in a restricted parameter space. The robustness of T— linear re- 
sistivity in the cuprates makes this type of scenario unlikely. What about new 
scenarios? Let us entertain the possibility that an additional length ^ is rele- 
vant and diverges as ^ oc with a > 1. In the calculation of the correlation 
volume in Eq. (jlj), one must replace with ^'^ ^ i'^ = Consequently, 
^(^/O) with h{y) = y^^ a general scaling function. In essence, one is reduc- 
ing the effective dimensionality such that d ^ d* = d — A(a — 1). T— linear 
resistivity obtains ii z = 2 — d*. The reduction in the effective dimensionality, 
A (a — 1), can now be flne-tuned so that (i* < 1, thereby resulting in physically 
permissible values of the dynamical exponent, z > 1. Such flne scripting of 
two length scales is also without basis at this time. 



3 Spectral Weight Transfer 

It is worth pointing out that non-commutative [24] fleld theories do permit 
z < 0. Central to z < is UV-IR mixing and as a consequence a break- 
down of the standard Wilsonian renormalization scheme. While a fleld theory 
based on non-commuting coordinates is undoubtedly not applicable to the 
cuprates, UV-IR mixing is certainly present as it has been documented ex- 
perimentally both in the normal state [25|26ll27f28f29j as well as the supercon- 
ducting state [30f3T][32] . Hence, the relevant question is: Can spectral weight 
transfer give rise to an additional length scale? 

To answer this question, we recount the well known argument [33|3l] on doping- 
dependent spectral-weight transfer across the Mott gap. Consider the half- 
flUed one-dimensional chain of one-electron atoms shown in Fig. ([2]) described 
by the Hubbard Hamiltonian, 
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Fig. 2. Spectral weight transfer in a doped Mott insulator. The Photoelectron spec- 
trum (PES) denotes the electron removal states while the electron-addition states 
are located in the inverse photo-electron spectrum (IPES). The on-site charging 
energy is U. Removal of a single electron results in the creation of two single par- 
ticle states at the top of the lower Hubbard band. By state conservation, one state 
comes from the lower and the other from the upper Hubbard band and hence spec- 
tral weight transfer across the Mott gap. 

i,j,cr i i(T 

in which electrons hopping on a lattice between neighbouring sites with am- 
plitude tij = taij and chemical potential /i pay an energy cost U anytime they 
doubly occupy the same site. The operator Ci„ (c|^) annihilates (creates) an 
electron on site i with spin a and rii^ is the occupancy on site i with spin a. 
When the hopping term vanishes, we can treat the half-filled system as having 
one electron per site. For a chain containing N sites, there are N ways to re- 
move an electron and ways to add an electron with an energy cost U. The 
operators that describe such excitations are = Qo-(1 — nj_o-) which adds an 
electron to an unoccupied site and rjj^ = cl^rii^^ which creates a doubly occu- 
pied site. Such operators create excitations in the lower and upper Hubbard 
bands, respectively, noted also in Fig. ([2]) as the photo-electron and inverse 
photo-electron spectra. When a single hole is created, the number of ways to 
remove an electron is now — 1. Surprisingly, the number of ways of adding 
an electron to the system so that the energy cost is U also decreases to — 1. 
This means that there are two less states at high and low energies. Where are 
the remaining two states? They correspond to the two ways of occupying the 
empty site with spin up and spin down electrons. Such states correspond to 
the addition part of the low-energy spectral weight and hence lie immediately 
above the chemical potential as illustrated in Fig. ([2]). One of the states must 
come from the upper Hubbard band (UHB) as the high energy part now has a 
spectral weight of iV — 1, the other from the lower Hubbard band. Hence, for 
a single hole, there is a net transfer of one state from high to low energy. In 
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general, simple state counting yields 2x for the growth of addition part of the 
low-energy spectral weight, A(x) and 1—x for the depletion of the high energy 
sector. In a Fermi liquid, adding holes simply creates quasiparticles near the 
chemical potential and hence cannot involve the high energy part. In actu- 
ality, the dynamical contribution [331IM] to the LESW results in A(x) > 2x. 
The dynamical LESW corresponds to virtual excitations to the UHB. Hence, 
in a strongly correlated system, some of the low-energy degrees of freedom 
arise from the high energy scale. A bonus that can be explained from spectral 
weight transfer is the sign change of the Hall coefficient seen widely [5^361137] 
in the cuprates. Particle-hole symmetry is restored in the lower Hubbard band 
when the addition and removal parts of the low-energy spectral weights are 
equal. Equating the atomic {U = oo) values for these weights, namely 1 — x 
for the addition spectrum and 2x for the removal part, we obtain that the 
maximum doping level at which the Hall coefficient changes sign is x = 1/3. 
This upper bound correlates with the experimentally observed value for the 
vanishing of the Hall coefficient [35|36|37j in the cuprates and is precisely the 
value obtained [38|l39] from strong-coupling calculations at large f/ ~ cxd. 

Spectral weight transfer described above can be reformulated in terms of the 
composite operators for the lower and upper Hubbard bands. Write the elec- 
tron operator as Qo- = + rjif^ and the corresponding spectral function. 



explicitly the imaginary part of the single particle Green function, contains 
diagonal terms as well as a mixing term which carries the information regard- 
ing the interconnectedness between the high and low energy scales. It is from 
the mixing term that spectral weight transfer arises. Shown in Fig. ([3]) is the 
mixing term computed using a self-consistent cluster method|40j). Regardless 
of filling, the cross term has both negative and positive contributions. This 
structure arises necessarily because the integral of over all frequency 



yields the equal time correlator {{C,ia,Via}) = 0; whose vanishing maintains 
the Pauli principle. This implies that is either zero, which it is not, or 
it must have both positive and negative parts, representing constructive and 
destructive interference, respectively, between different regions in energy. At 
half-filling, upper panel, if particle- hole symmetry is present, the contribu- 
tions of the cross term below and above the chemical potential sum to zero 
independently, indicating that all energy scales need not be retained to ensure 
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Fig. 3. Cross correlation or quantum interference A^^ between the upper and lower 
Hubbard bands at half-filling, n = 1 and at n = .95 at T = O.li. The dip at the 
chemical potential in the lower panel represents the pseudogap. The inset shows 
that this dip leads to a vanishing density of states at zero temperature and hence 
an orthogonality catastrophe. 

the Pauli principle. Such is not the case, however, at finite doping. The lower 
panel in Fig. ([3]) indicates that a pseudogap develops at the chemical poten- 
tial, indicating an orthogonality catastrophe. Because a pseudogap subtracts 
spectral weight at low energies and transfers it to intermediate to high energy, 
the sum rule which ensures the Pauli principle is satisfied only when is 
integrated over all energy scales not simply up to the chemical potential (or 
some intermediate energy cutoff) as would be the case in projected models. 

That projection changes fundamentally the statistics of the particles arises 
from the truncation of the Hilbert space. Consider the standard way of formu- 
lating the t — J model from the Hubbard model: block diagonalise the Hubbard 
model into sectors with a fixed number of doubly occupied sites via a similar- 
ity transformation 131], S, which connects sectors that differ by a single double 
occupancy. In such a scheme, the electron operator in the low-energy sector is 
transformed to 



where P removes double occupancy. To lowest order in t/U, ata = and 
^h^ii ~ 0- That is, projected operators in Eq. (fTTl) block two states from being 



9 




<. « — ^ 



() e e p- 



T/t 



Fig. 4. Resistivity as a function of temperature for the Hubbard model (with 
U = lot) using the spectral function computed previously by Stanescu and 
PhillipsgT] for fillings n = 0.97,0.95,0.9,0.85,0.80. po = h/e"^. The inset shows 
the doping dependence of the pseudogap energy scale, T*. The dependence obeys 
the functional form, J(l — x^^Q/4a^), where ^do is the average distance between 
doubly occupied sites, the percentage of which is indicated in the table. 

occupied rather one as would be the case for fermions. Simply, projection of the 
electron onto the lower Hubbard band does not result in a fermion. In fact, such 
projected objects are exclusons obeying generalized exclusion statistics |42]. 

Since truncating the Hilbert space from four to three states per site changes the 
statistics, the limits U ^ oo and L oo do not commute. Whether or not this 
changes the physics fundamentally is open for debate. However, if transport 
is governed by a finite length scale for double occupancy, ^do, then the correct 
order of limits is L — > cxd followed by f/ — > oo. Projected models such as the 
t — J model correspond to the opposite regime, f/ — *• cx3 and then L — > cxo. 
In this case, ,^do = oo. K finite value of ^do would certainly function as the 
additional length scale that is required to describe T-linear resistivity. In recent 
cluster calculations [S] on the Hubbard model, we observed an insulating state 
(see Fig. (jlj)) below a critical doping level as is seen experiment ally p1l3]|^ll5]|6] . 




The pseudogap arises [H] from the eneryg gap between local cluster states 
with total spin S and S — 1. This energy difference is essentially the triplet- 
singlet splitting and hence scales as /U . We argued that the insulating state 
persists as long as rihi'io < -^^5 = xi^L/af the number of holes, rihi'io = -^^ 
defines the percolation limit. By calculating the percentage of doubly occupied 
sites, we obtained ,^do numerically and plotted the T*-line, J(l — cxi^^^o/ a)"^), 
in the inset of Fig. (jl]). The agreement of this phenomeno logical fit with the 
resistivity data in which a metallic state obtains at a; = 0.1 lends credence to 
our assertion that ^do is the relevant length scale for the pseudogap. Similar 
calculations on the t — J model find a metallic state [ISHG] consistent with the 
lack of commutativity of the limits U ^ 00 and L — > 00. But this question is 
far from settled. 
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4 Unifying Principle: Zeros 



To help resolve the potential non-commutativity of t/ — oo and L — > oo, 
we seek a unifying principle from which spectral weight transfer, pseudogap 
phenomena, and broad spectral features emerge. Such a principle does exist 
and stems from the analytical structure of the single-particle Green function 
whenever a gap is dynamically generated. In Fermi liquids, the single-particle 
Green function has poles at the quasiparticle excitations. The divergence of the 
single-particle Green function defines the Fermi surface. In insulators, no such 
divergence obtains because quasiparticles are absent. However, Luttinger's|48] 
theorem 



has been argued [19] to be still valid even when quasiparticles are absent. The 
essence is that the integral defining the particle density is performed over 
a slice in momentum space where the Green function is positive. While nu- 
merous approximate calculations [51] have found surfaces of zeros in ID [51] or 
quasi- ID [52|53] systems, it is unclear if the zeros are an artifact of the approx- 
imations. Further, it is implicitly assumed in the work of Dzyaloshinskii[4^ 
that the volume of the surface of zeros equals the particle density. Should 
this not be the case, Eq. f[T^ must be appended with a correction to recover 
the particle density. In fact, this asumption has been questioned for a specific 
model for a Kondo insulator |50j. 

Recently, we have found a general proof [5lj that allows us to identify precisely 
where the surface of zeros occurs regardless of spatial dimension provided 
certain symmetries obtain. We expand on that proof here. To proceed, we 
note that the causal nature of the retarded Green function 



permits it to be constructed entirely from the spectral function (that is, the 
imaginary part) 



through the standard Hilbert representation. For a Mott insulator, a gap of 
order U occurs in the spectral function. We will take the gap to have a width 
2A centered about 0. Within the gap, A(k, e) = 0. This is a necessary condition 
for any gap. Consequently, in the presence of a gap, the real part of the Green 
function evaluated at the Fermi energy reduces to 





(13) 
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We now prove that when particle- hole symmetry is present, the retarded Green 
function is an even function of frequency at the non-interacting Fermi surface. 
As a consequence, Eq. f fTSj) is identically zero along that momentum surface. 
To proceed, we consider the moments 



M^ik)^ l^u-duGTiku) (16) 

of the Green function. For simplicity, we have set h = 1. Using the Heisenberg 
equations of motion, we reduce the moments in real space 



1 r 



J n times; 



+ {{c,^,[---[cl,H]---H],HU^,,}) 



(17) 



to a string of commutators of the electron creation or annihilation operators 
with the Hubbard Hamiltonian. The right-hand side of this expression is eval- 
uated at equal times. To evaluate the string of commutators, it suffices to 
focus on the properties of 



ir£) = [---[Q.,i/],---i/]ntimes, (18) 

where by construction, K^^^ = Cia- We write the Hubbard Hamiltonian as 
H = Ht + Hu where Hu includes the interaction as well as the chemical 
potential terms. The form of the first commutator. 



3 

suggests that we seek a solution of the form 



^r = Ew. (20) 

i 

where Q-"^ = [•■• [cjo-, -f^;/], ■■■ i^;7]n times involves a string containing Hjj n 
times and in Aj^, Ht appears at least once. Our proof hinges on the form 
of qI"'* which we write in general as Q'f^ = anCicjUi^^ + (3nCia- The solution 
for the coefficients 
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= (t/-/i)a„ + ?7(-/x)" 

[Zi) 

is determined from the recursion relationship Qjv^^'' = [Qi^\ Hu]- In the mo- 
ments, the quantity which appears is 



{{Qt\ J) = 5ij + ^ 5,,-7„. (22) 

Consequently, the moments simplify to 

M^{^,J) = 5,,7„ + lT.t^l c]J) + h.c.) . (23) 



The criterion for the zeros of the Green function now reduces to a condition on 
the parity of the right-hand side of Eq. fl23l) . Consider the case of half-filling 
and nearest-neighbour hopping. Under these conditions, (njo-) = 1/2 and by 
particle- hole symmetry, /i = U/2. The expressions for an and Pn lays plain 
that the resultant coefficients 



..-{^r-±^ (24, 

vanish for n odd. Consequently, Ga{k,u) is an even function if the second 
term in Eq. fl23l) vanishes. In Fourier space, the second term is proportional 
to the non-interacting band structure, which in the case of nearest-neighbour 
hopping is t{k) = — 2tY,i=i cos ki. ^ result, the Green function only has 
even moments at the momenta for which X^iLi cos fcj = 0. The vanishing of 
t{k) defines the non-interacting Fermi surface. Consequently, the surface of 
zeros is pinned at the Fermi surface of the non-interacting system at half-filling 
whenever a Mott gap opens in the presence of particle-hole symmetry. This 
constitutes one of the few exact results for Mott insulators that is independent 
of spatial dimension. The only condition for the applicability of our proof is 
that the form of the gap leads to finite integrals in Eq. f|T5l) . Hence, the minimal 
condition is that the density of states vanishes at zero frequency as w", a > 0. 
Note, this encompasses even the Luttinger liquid case in which < a < 1. 

To prove that this observation applies regardless of the range of the Coulomb 
interactions, we consider an argument directly from particle-hole symmetry. 
Under a general particle-hole transformation. 



e*Q-'4, (25) 
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with the chemical potential fixed at /i = ?7/2 in the Hubbard model, the 
spectral function becomes 



A„(k,uj) = A^{-k- Q + 2mi,-uj), (26) 

Hence, the spectral function is an even function of frequency for k = Q/2+n7r. 
Consider one dimension and nearest-neighbour hopping. In this case, the sym- 
metry points are ±7r/2, the Fermi points for the half-filled non- interacting 
band. In two dimensions, this proof is sufficient to establish the existence of 
only two points, not a surface of zeros. To determine the surface, we take 
advantage of an added symmetry in higher dimensions. For example, in two 
dimensions, we can interchange the canonical x and y axes leaving the Hamil- 
tonian unchanged. This invariance allows us to interchange and ky on the 
left-hand side of Eq. (!26l) resulting in the conditions 

ky = —kx — q + 2n7r (27) 
and by refiection symmetry 



— ky = —kx — q + 2n7r, (28) 

where Q = (g, g). For nearest neighbour hopping, the resultant condition, 
kx ^ ky = — TT + 2?7,7r, is the solution to cos kx + cos ky = 0, which defines 
the Fermi surface for the non-interacting system. When nearest-neighbour 
Coulomb interactions are present, the chemical potential changes to {U+2V)/2 
at the particle-hole symmetric point and the proof follows as before. Arbitrary 
density-density interactions and hence all Coulomb interactions are indepen- 
dent of Q under a particle-hole transformation. Hence, our proof applies re- 
gardless of the range of the Coulomb interactions. A proposed proof[56j for the 
existence of zeros in ID relies explicitly on rotational invariance of the Green 
function. Such an invariance is not applicable to the Hubbard model as the 
charge and spin velocities are necessarily different as long asU ^ 0. Hence, our 
proof constitutes the only proof for zeros exist in the Hubbard model regard- 
less of spatial dimension. Finally, existing calculations based on the random 
phase approximation in which zeros were found at the non-interacting Fermi 
surface in ID [51] and quasi- 1 D [52][53] systems are consistent with our proof. 

While this proof applies strictly at half-filing with particle-hole symmetry, it 
places constraints on how the Fermi surface can evolve at finite doping in 
Mott systems. The consequences of this result follows. First, in the 

presence of the Mott gap, the self-energy diverges at zero frequency along the 
surface of zeros. To prove this, observe that the single-particle Green function 
can be written as Ga{k, u) = — e{k) — 'RT.aik, u) — i'^'E^j^k, a;)), where 
S is the self-energy. Near the Luttinger surface, G oc k — kj^. Q'S(/c,0) oc 
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S{k — ki). Note, however, that QG{k,uj) = for all energies within the gap. 
Now simply invert the Green function to obtain that 3fJEo-(fc, uj) (x {k — ^l)"^, 
thereby proving the assertion. Such a divergence is a clear indication that 
perturbation theory breaks down. Zeros are a concrete manifestation of this 
breakdown. The divergence of 3fJS prevents the renormalized energy band, 
e{k) + 3fJScr(/c, 0) from crossing the chemical potential, thereby resulting in an 
insulating state. This reinterpretation of the Mott insulating state provides 
a general way of understand how insulating states arise through dynamical 
generation of a gap in the absence of symmetry breaking. In dynamical cluster 
calculations, the divergence of S(7r, 0) [S7j has been observed but attributed 
to antiferromagnetism. What our proof makes clear is that the divergence of 
the self-energy is a result of Mottness itself; that is, it is not concomitant to 
T = ordering. 

Second, the volume of the surface of zeros equals the particle density only 
when particle-hole symmetry is present. In proving this, we recall the most 
general way of writing the electron density 



7 = --. + 2y^^G(k,.)^-^. (29) 

Under Eq. ([25]), G{k,uj) -^(-k - Q, -cu) andE(k,w) ^ -S(k-Q,-cu), 
our Hubbard Hamiltonian remains invariant with the chemical potential fixed 
at the particle-hole symmetric value of /i = U/2, and the additional contribu- 
tion to the electron density in Eq. flSU]) vanishes identically. For Fermi liquids, 
in which the self-energy is smooth, the second term in Eq. flSU]) is integrated 
by parts [?] 



Jd'^kjdu (G(k,.)^^ 
^ 

-G(-k-Q..) ^^'-';-Q-"' )=0 (30) 

and shown to be identically zero. For a Mott system lacking particle-hole 
symmetry but possessing a divergent self-energy, the standard integration by 
parts [?] of the second term in Eq. fl30l) is not valid. The singular part of Gd^jT, 
will always integrate to a non-zero value. That is, in the absence of particle- 
hole symmetry, no sum rule allows us to equate the volume of the surface of 
zeros with the particle density as argued for perturbatively in a specific model 
for a Kondo insulator |50j. 

Third, Zeros are absent from projected models at half-filling. Under projec- 
tion, the real part of the Green function reduces to the first integral in Eq. 
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( ITSi) because projection does not preserve the contribution above the gap. 
This integral is always positive and hence zeros are absent. Transforming the 
operators in the t — J model to respect [58j the no double occupancy condi- 
tion is of no help as the problem stems from the loss of spectral weight above 
the gap once projection occurs. This result points to asymptotic slavery [lO] 
of the Hubbard model or equivalently an example of the non-commutativity 
discussed earlier in that a low-energy reduction changes the physics. Stated 
another way, the t — J model violates Luttinger's theorem at half-filling. The 
essence of this breakdown is that there is no kinetic energy in the t — J model 
at half-filling, simply interacting spins. If there is no kinetic energy, there can 
be no surface of zeros. Such is not the case in the Hubbard model since no 
eigenstate of the Hubbard model has a fixed number of doubly occupied sites. 

Fourth, even at infinitesimal doping, the t-J and Hubbard models are not 
equivalent. Strictly speaking, at half-filling, it is not entirely appropriate to 
compare the t — J and Hubbard models because the former has no spectral 
weight above the chemical potential. Perhaps the proper way is to compare 
the two models in the limit of infinitesimal doping (see Figs. ([5^) and (c)). 
For one hole (n = numerical and analytical studies [5^6U|61] find a 
quasiparticle in the t-J model with weight J/t at (7r/2,7r/2) whereas in the 
Hubbard model[62], the quasiparticle weight vanishes as Z oc 6* > 0, L 
the system size. This result is consistent because the one- hole system should 
closely mirror the half-filled system, which for the Hubbard model must have 
a surface of zeros for n = 1 while no such constraint exists in the t-J model. 
This has profound consequences for the minimal model needed to describe the 
cuprates. 

Fifth, band structure cannot affect the existence of zeros in the presence of 
the Mott gap. At present, our proof apphes to any kind of band structure 
that is generated from hopping processes which remain unchanged after the 
application of Eq. fl2^ . In general, the two kinds of hopping processes trans- 
form as e(7r — k^, n — ky) = —e{kx, ky) and t'{n — /c^;, vr — ky) = t'{kx, ky). The 
latter is relevant to the cuprates. If only such hopping is present, the surface of 
zeros is no longer the diagonal (vr, 0) to (0, vr) (or the point ±7r/2 in ID), but 
rather the "cross" (0,7r/2) to (7r,7r/2) and (vr/2,0) to (7r/2,7r) (or, in ID, the 
points ±7r/4 and ±37r/4). When both types of hopping are present, no sym- 
metry arguments can be made. However, at this point it is imperative that 
we distinguish between the ID and 2D cases. In ID, if only two Fermi points 
exist, the Luttinger theorem requires that the zeros be at the symmetric Fermi 
points, ±7r/2. But in principle there can be 2n Fermi points in ID. Consider 
the case of t and t' hopping. Four Fermi points exist for t' > t/2. In such 
cases, the zeros need not be at the Fermi surface of the non-interacting sys- 
tem. Such ID cases share the complexity of the general case of an asymmetric 
band structure in 2D. When this state of affairs obtains, we can establish the 
existence of zeros by a key assumption: the Green function is a continuous 
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function of the hopping parameters t and t'. When only t is present, Ra{k,0) 
has one sign (plus) near k = (0, 0) (or, in ID, k=0), and the opposite (minus) 
near k = (vr, vr) {k = it in ID) and will vanish on the zero line. Alternatively, 
if we have t' hopping, R„{k,0) will have a certain sign near k = (0,0) and 
k = {re, it), and the opposite sign near k = (0,7r) and k = (tt, 0) and will 
vanish on the "cross". From continuity, for t' <^ t, Rcr{k]t,t') will have the 
same sign structure as Ra{k; t, t' = 0). That is, it will change sign when going 
from (0, 0) to (vr, vr) regardless of the path taken. Therefore, the line of zeros 
exists for small enough t', the relevant limit for the cuprates. In the opposite 
limit, t' ^ t, a. similar argument holds. 



Finally, the surface of zeros defines the pseudogap at finite non-zero dop- 
ing. At finite hole-doping [x = 1 — n), the chemical potential jumps into the 
lower Hubbard band. For a; ~ 1, R(j{k,0) must vanish on some momentum 
surface because most of the spectral weight at (vr, vr) still resides above the 
chemical potential, whereas at (0, 0) it resides below. Consequently, a surface 
of zeros should still be present so long as a gap, commonly known as the 
pseudogap, opens up to some doping Xc- Satisfying the zero condition, Eq. 
(fT5l) . requires spectral weight to lie immediately above the chemical potential. 
Spectral weight transfer across the Mott gap depicted in Fig. ([2]) and dis- 
cussed extensively previously mediates the zeros. The pseudogap that results 
is necessarily asymmetrical for small x because the 2x peak must lie closer 
to the chemical potential than the 1 — a; contribution in order for Eq. f|T5l) 
to vanish. Hence, an asymmetric pseudogap is a direct indication that Mot- 
tness rather than weak-coupling/ordering physics is the operative cause. In 
fact, such an asymmetry has been observed in the pseudogap regime in the 
cuprates [55]. The real question is how do the zeros evolve upon doping. As 
shown in Fig. (Q, either the surface of zeros vanishes abruptly above a critical 
doping or it shrinks smoothly. In either case, the pseudogap must vanish at 
a finite number of fc— points. A hard gap would involve all momenta. Hence, 
the pseudogap must vanish as tu", with a > . Because poles cannot exist 
within an e neighbourhood of a zero if the pseudogap arises from Mottness, 
quasiparticles cannot survive at the nodes. Weak-coupling or Fermi liquid- 
type scenarios permit nodal quasiparticles. Numerics |64|65ll66j point to such 
a non-symmetry breaking pseudogap in the Hubbard model at finite doping 
whose cause is Mottness. Also, Stanescu and Kotliar[n7| have shown con- 
vincingly using a cellular dynamical cluster method on the Hubbard model 
with nearest-neighbour hopping that the spectral features of a momentum- 
dependent pseudogap are directly related to the surface of zeros. This is the 
first extensive numerical study that has established the hard connection be- 
tween zeros and the pseudogap. 
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Fig. 5. Evolution of the surface of zeros in the first quadrant of the FBZ. Yellow 
indicates KG > while blue < 0. The Hubbard model is constrained to have a 
surface of zeros as n — > 1 whereas the t-J model is not. The two options upon doping 
represent weak-coupling (la or c to d) and strong-coupling (la to b). The transition 
from (a) to (b) requires a critical point at ric whereas (d) does not. Experiments [2|6] 
indicating an insulating state for n > ric are consistent with an abrupt transition 
from (a) to (b). 



5 Concluding Remarks 



Because of its explanatory power, the surface of zeros can be argued to play a 
role in doped Mott insulators analogous to the Fermi surface in Fermi liquids. 
However, the only feature they share in common, is their volumes when the 
Hamiltonian has particle-hole symmetry. In Fermi liquids, no spectral weight 
transfer occurs, no pseudogap phenomena exists and quasiparticles are well 
defined. Just the opposite obtains in doped Mott systems. That broad spectral 
features naturally follows from the surface of zeros is clear because sharp quasi- 
particles require the renormalized energy band to cross the chemical potential. 
This is not possible as long as ^J2a{kL, 0) diverges as has been demonstrated. 
Hence, one of the most nettling problems with the cuprates, namely broad 
spectral features [8] has a natural resolution on the surface of zeros. The ze- 
ros also offer a concrete realisation of the non-commutativity of ?7 — >■ cxo and 
L — > cxD as they are present in the Hubbard model at half-filling but not in 
the projected t—J model. A vanishing of the surface of zeros above a critical 
doping level will define a length scale associated with the gap. Such a length 
scale will undoubtedly play a role in the transport properties in the strange 
metallic regime. A complete treatment of the transport properties and the 
nature of the excitations on the surface of zeros will necessitate a complete 
field theory of Mottness. 
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